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Computational and experimental analyses of a PICS — Pilot-In-Can-Swirler technology 
injector, developed by United Technologies Research Center (UTRC) are presented. NASA 
has defined technology targets for near term (called “N+l”, circa 2015), midterm (“N+2”, 
circa 2020) and far term (“N+3”, circa 2030) that specify realistic emissions and fuel 
efficiency goals for commercial aircraft. This injector has potential for application in an 
engine to meet the Pratt & Whitney N+3 supersonic cycle goals, or the subsonic N+2 engine 
cycle goals. Experimental methods were employed to investigate supersonic cruise points as 
well as select points of the subsonic cycle engine; cruise, approach, and idle with a slightly 
elevated inlet pressure. 

Experiments at NASA employed gas analysis and a suite of laser-based measurement 
techniques to characterize the combustor flow downstream from the PICS dump plane. 
Optical diagnostics employed for this work included Planar Laser-Induced Fluorescence of 
fuel for injector spray pattern and Spontaneous Raman Spectroscopy for relative species 
concentration of fuel and C0 2 . 

The work reported here used unheated (liquid) Jet- A fuel for all fuel circuits and cycle 
conditions. The initial tests performed by UTRC used vaporized Jet-A to simulate the 
expected supersonic cruise condition, which anticipated using fuel as a heat sink. 

Using the National Combustion Code a PICS -based combustor was modeled with liquid 
fuel at the supersonic cruise condition. All CFD models used a cubic non-linear k-epsilon 
turbulence wall functions model, and a semi-detailed Jet-A kinetic mechanism based on a 
surrogate fuel mixture. Two initial spray droplet size distribution and spray cone conditions 
were used: 1) an initial condition (Lefebvre) with an assumed Rosin-Rammler distribution, 
and 7 degree Solid Spray Cone; and 2) the Boundary Layer Stripping (BLS) primary 
atomization model giving the spray size distribution and directional properties. Contour and 
line plots are shown in comparison with experimental data (where this data is available) for 
flow velocities, fuel, and temperature distribution. The CFD results are consistent with 
experimental observations for fuel distribution and vaporization. 

Analysis of gas sample results, using a previously -developed NASA NOx correlation, 
indicates that for sea-level takeoff, the PICS configuration is predicted to deliver an EINOx 
value of about 3 for the targeted supersonic aircraft. Emissions results at supersonic cruise 
conditions show potential for meeting the NASA goals with liquid fuel. 
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Nomenclature 


AST 

= 

Advanced Subsonic Technology 

CFD 

= 

computational fluid dynamics 

CFL 

= 

Courant-Friedrichs -Le wy condition 

c P 

= 

specific heat 

CPU 

= 

central processing unit 

8 

= 

turbulent dissipation 

EINOX 

= 

emission index for oxides of nitrogen 

ERA 

= 

environmentally responsible aviation 

FAR 

= 

fuel-to-air ratio 

FWHM 

= 

full width at half maximum 

JST 

= 

Jameson-Schmidt-Turkel dissipation scheme 

k 

= 

turbulent kinetic energy 

LTO 

= 

landing-takeoff cycle 

V 

= 

molecular vibrational energy 

NCC 

= 

National Combustion Code 

PICS 

= 

Pilot-In-Can- Swirler 

P3 

= 

combustor inlet pressure 

PLIF 

= 

planar laser-induced fluorescence 

PLS 

= 

planar laser scatter 

T3 

= 

combustor inlet temperature 

SLTO 

= 

sea level takeoff 

SRS 

= 

spontaneous Raman spectroscopy 

Vdrop 

= 

liquid phase droplet group velocity vector 

Vgas 

= 

gas phase velocity vector 

V' 

Y gas 

= 

turbulent fluctuations of the gas phase velocity vector 


I. Introduction 

F OR more than 40 years, NASA has sustained programs to reduce the environmental effects of aviation. A 
major focus of these programs has been reducing the emissions of nitrogen oxides (NOx). NOx emissions 
decrease the protective ozone layer in the stratosphere and increase smog and ozone in the lower troposphere 1 . To 
prevent damage to the protective ozone layer, NASA programs have focused on reducing NOx emissions at cruise 
conditions for supersonic flight. To reduce the emissions of NOx in the lower troposphere, NASA programs have 
also focused on reducing NOx emissions during the landing -takeoff cycle in subsonic flight. 

In addition to decreasing NOx emissions, NASA has also addressed the reduction of carbon dioxide emissions 
by increasing engine fuel efficiency. Increased fuel efficiency is typically achieved by increasing the engine 
operating pressure ratio, which increases combustor inlet temperature; however, NOx formation rates increase with 
higher temperatures. Without an improvement in combustor technology, higher efficiency engines will have higher, 
not lower, NOx emissions: improved low-NOx combustor technologies need to be developed. NASA has defined 
technology targets for the near term (called “N+l”, circa 2015), midterm (“N+2”, circa 2020) and far term (“N+3”, 
circa 2030) that specify realistic emissions and fuel efficiency goals for commercial aircraft. In the far term, the 
emissions goal for supersonic cruise is less than 5 g/kg of fuel burned, with fuel efficiency targeted at between 3.5 - 
4.5 passenger-miles per pound of fuel 2 . The Supersonics Project of the NASA Fundamental Aeronautics Program is 
charged with managing the challenge to develop the technologies required to meet these goals. The technical 
approach taken incorporates using physics -based modeling for combustion CFD to improve supersonic cruise 
emissions predictions, and using validation experiments to assure the models work. These models are applied to 
promising low emission concepts to aid in design and development. NASA has partnered with engine and fuel 
injector manufacturers to develop practical technology to meet the stringent emissions goals requirements. 

One such partnership is with United Technologies Research Center (UTRC), which devised a fuel 
injection/mixing concept to meet the far term criteria using a fuel-lean approach compatible with the cycle 
conditions devised for a Pratt & Whitney supersonic transport notional engine. The concept is called the Pilot -In- 
Can S wirier (PICS). In the PICS injector concept, each swirler contains a pilot “can” concentrically located inside 
the main stage swirler (which dumps into an annulus), and in which the pilot stage flame is isolated from the main 
stage flame to help reduce interaction between zones for improved turndown capability. One key aspect of the 
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concept is that the fuel for the main fuel stage is injected as a gas 3 . Gas phase fuel air mixing theoretically has the 
benefit of reducing the additional mixing time that would have been required for fuel atomization and vaporization. 
Additionally, if the fuel is able to be preheated by using it as a heat sink for cooling the aircraft, other efficiencies 
can be gained. 

Previous evaluation of the PICS injector performed at UTRC 3 included comparison of cold-flow gas-gas mixing 
measurements with analytical predictions for supersonic cruise conditions, and combustion tests at both supersonic 
cruise and idle conditions (pilot only) using vaporized Jet- A fuel for the main stage and non-heated, liquid Jet- A for 
the pilot. The gas-gas mixing results showed acceptable variation in mixture fraction at the primary swirler exit. The 
emissions results at supersonic cruise satisfied the NOx emissions target with combustion efficiency > 99.99%, and 
at idle conditions also had high combustion efficiency with EINOx < 6. 

This paper reports results from CFD analysis and combustion testing conducted at the NASA Glenn Research 
center, to explore PICS injector performance when using liquid (not gaseous) fuel. Testing of the PICS injector at 
NASA used unheated liquid fuel for all fuel circuits and all power settings, and expanded the previous test matrix. 
For these liquid-fuel tests at NASA, UTRC designed and fabricated a PICS injector having a reduced number of 
main fuel -injection orifices as compared to the earlier vaporized-fuel design. The smaller number of injection 
orifices ensured adequate fuel -injection velocity when using liquid fuel, which has a lower volume flow rate than 
vaporized fuel because liquid is more dense than vapor. Optical diagnostics techniques were applied at the 
combustor inlet to examine injector performance for fuel liquid versus vapor distribution and mixedness, and 
exhaust gas analysis was applied at the combustor exit to measure emissions. The test results are compared with 
CFD results produced using the National Combustor Code (NCC), specifically at supersonic cruise conditions. 


II. PICS Hardware 

The PICS injector concept, illustrated in Fig. 1, consists of a pilot stage and a main stage, with turndown 
controlled by shifting the fuel distribution between the stages. To achieve this, the injector contains a pilot ’’can” 
centrally located within the main swirler, in which the pilot flame is isolated from the main-stage flame. When 
installed within a single-annular combustor dome, each injector is identical and receives the same fuel flow 
allocation. All tests and computations at NASA used a single PICS swirler, installed in a single 4 -inch by 4-inch 
sector. The single PICS injector hardware is shown in Fig. 2 and the computational envelope used for CFD is shown 
in Fig.3. 



Figure 1. PICS swirler-injector Figure 2. Photograph of the single 
concept. PICS injector, from aft-looking- 


forward view. 



Figure 3. Computation envelope 
used by the National Ccombustion 
Code for CFD simulations of PICS. 
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III. Experiment Setup 


A. PICS Installation 

The combustor subcomponent facility at NASA Glenn supplies heated, non-vitiated air to the test rig. The test 
rig uses a castable ceramic to form the combustor “liner”, so the four inch by four inch PICS module had adiabatic 
walls, rather than a water-cooled liner as used for the UTRC tests. Another difference in the NASA installation was 
that the test stand on which the PICS injector was installed is best suited for higher air and fuel flow rates. Because 
the two facilities are designed to accommodate different flows, the two sets of tests could not completely match the 
Pratt & Whitney supersonic cycle conditions explored in this study. Table 1 shows the Pratt & Whitney supersonic 
N+3 cycle conditions. Table 2 those lists the conditions used at NASA Glenn that incorporate the laser diagnostics 
measurements. 


Table 1. Pratt & Whitney supersonic N+3 cycle. 3 



P3, psi 

H 

u> 

o 

FAR/FARslto 

Supersonic cruise 50,000 ft, Ml. 8 

174 

1087 

1.10 

SLTO 

329 

890 

1.00 

Climb-66% 

235 

767 

0.82 

Approach-32% 

149 

634 

0.66 

idle 

78 

475 

0.58 


Table 2. Inlet condtions used for PICS tests with laser diagnostics. 


Nominal cycle 

P3, psi 

H 

u> 

o 

FAR/FARslto 

Supersonic cruise, N+3 

174 

975 

1.24 

~ Subsonic cruise, N+2 

250 

1000 

0.85 

Approach, subsonic, N+2 

205 

716 

0.75 

~ Idle, subsonic, N+2 

100 

425 

0.85 


B. Laser and Optical Diagnostics 

To determine the fuel spray pattern, we acquired 
two-dimensional images of fuel via planar laser- 
induced fluorescence (PLIF) and planar laser scatter 
(PLS); to measure major combustion species (N 2 , 
0 2 , CO, C0 2 , fuel, and H 2 0), we acquired one- 
dimensional images using Spontaneous Raman 
spectroscopy (SRS). We also acquired 
chemiluminescence images of Swan band C 2 ., 
which are visible in the reaction zone and show 
where carbon-carbon electronic transitions occur. 
(CH and C 2 excited species emit light in the visible 
region of a hydrocarbon flame spectrum.) Setup and 
data acquisition are described next. 



Figure 4. Schematic drawing that shows the key 
components used for planar laser-induced fluorescence 
and planer laser scattering. 
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The fuel PLIF, PLS, and C 2 * images were obtained using the same receiving optics and intensified CCD (ICCD) 
camera. A second set of fuel PLIF images was gathered using a second ICCD camera, located on the opposite side 
of the combustor. The setup is similar to that illustrated in Fig. 4, in which the cameras were positioned on the sides, 
to image light from an angle perpendicular to the combustor flow direction. A remotely controlled filter wheel 
allowed us to select the appropriate filter (FWHM of 10-nm, typical) to pass light through a UV-grade, f = 105-mm, 
f/4.5, macro camera lens. The light was then focused onto a gated, 16-bit, lk x lk pixel array, ICCD camera having 
a Gen II Super-Blue-Slow-Gate intensifier. The gate time used was 100-ns. The laser used was a 10-Hz, frequency- 
doubled Nd:YAG-pumped dye laser/frequency mixer system set up to produce wavelengths around 282-nm. The 
laser beam was formed into a sheet using a pair of cylindrical lenses, to obtain a sheet approximately 300 pm thick. 
We typically acquired images with on-chip averages of 200 gates, and traversed across the flow along the Y axis in 
1-mm increments. The laser sheet and collection optics were traversed together so as to maintain focus on the laser 
sheet. For PLIF, we collected fluorescence using a filter centered near 334-nm on one camera and near 313-nm on 
the other. For the elastically- scattered light of PLS, we used a filter centered near 280-nm for each camera. For 
chemiluminescence, we did not use the laser and collected light when positioned at Y=0 using a filter centered at 
514-nm. We collected PLIF, PLS and C 2 data for all test matrix conditions. 

The spontaneous Raman species measurements were obtained using an optical arrangement similar to that shown 
in Fig. 5. In this case, 532-nm light from a frequency-doubled Nd:YAG laser operated at 15-Hz provided the light- 
scattering source. The probe volume was formed using a spherical lens, to provide a narrow, high -energy-density 
region from which molecules can scatter light (the SRS signal is -10-6 smaller than for laser induced fluorescence). 
The shape of the probe volume was defined by a combination of input and collection optics which produces a 
cylinder 6.67-mm high — aligned in the vertical (z-) — with approximately a 2-mm diameter. The scattered light from 
the molecules in the probe volume passed through a 550-nm long pass filter and was collected using an f = 60-mm 
spherical lens coupled to a vertically-oriented, linear fiber array, consisting of thirty-seven, 200-pm-diameter fibers. 
The fiber transmitted the light to an f/4, f = 300-mm spectrometer with a 600 groove/mm grating blazed at 500 nm. 
The light dispersed by the grating was focused onto 
a High Q Blue, gated, intensified camera having an 
array size 1024-pixel- wide by 256-pixel- high. The 
camera CCD pixel size combined with the 600 
line/mm grating produced a -0.1667 nm/pixel 
spectral resolution. The camera array was binned in 
the spatial (vertical) dimension into five rows, each 
40 pixels high, with the remaining pixels not used. 

These five rows provide for simultaneous 
measurement at five discrete locations in the z- 
direction. Sets of 100, 10-gate on-chip averages 
were collected at each position traversed in x and y. 

The five binned rows created five spectra in which 
the peaks were used to identify the relative species 
concentrations. The use of Raman was restricted to 
those inlet conditions expected to have little or no 
liquid fuel so that the interference from Mie scatter 
is minimized; Mie scatter has a signal roughly a 
million times greater than SRS. 

For all measurements, there may be a slight 
difference in the true spatial coordinates because the heated inlet air causes the test rig will to shift axially (x-), 
horizontally (y-) and rotationally about its centerline. We took measures during the tests to account for translations, 
but the rotation was hard to quantify and not addressed. 



Figure 5. Illustrations that shows the key components 
used for spontaneous Raman scattering experiments. 


C. PICS Modeling 

7. The National Combustion Code 

The National Combustion Code (NCC) is a state of the art CFD program specifically designed for combustion 
processes. A short summary of the features of the NCC pertaining to this paper are: the use of unstructured grids 4 , 
massively parallel computing - with almost perfectly linear scalability 5 ’ 6 on non-spray cases up to four thousand 
central processing units (CPU), a dynamic wall function with the effect of adverse pressure gradient 7 , low Reynolds 
number wall treatment 8 , and a cubic non-linear k-epsilon turbulence model 9, 10 , lagrangian liquid phase spray 
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model 11 , and stiff laminar chemistry integration. Recently, viscous low-speed preconditioning 12,13 has been added to 
improve the low-speed convergence of the NCC in viscous regions. The combination of these features is usually not 
available in other CFD codes and gives the NCC an advantage when computing recirculating, turbulent, reacting, 
spray flows. Previously, the NCC has undergone extensive validation studies for simple flows 14 , complex flows 15 , 
NO x emissions prediction performance 16 , and traditional gas turbine combustor/injectors 17 . 

2. Geometry and Mesh Generation 

UTRC created the orginial PICS solid model using NX 18 and transferred it to NASA in ParaSolid format as an 
“air-solid”, the solid solid geometry that only represented the fluid flowing through the PICS combustor and NASA 
test cell. The file was then imported into SolidWorks, “regenerated” and exported as an ACIS® file. The ACIS file 
was then imported into Cubit 19 for automatic geometry clean up and manual simplification, while also removing 
surface artifacts (like holes where tubes once existed) from the simplification process. 

The single, simplified “air-solid” was then imported as an IGES file into Pointwise 20 . Because ignition transients 
tended to greatly increase wall clock computational time as mesh size increased, an optimized mesh of 5.5 million 
isotropic tetrahedral elements was generated. We believe this mesh balanced accuracy and reasonable solution time 
for this engineering effort. Non-reacting simulations were run as a check using a stretched tetrahedral mesh of 18 
million elements. The 5.5 million element matched the global pressure drop to within 20% of the 18 million element 
mesh value. 

Chemistry Modeling 

Ideally, we would prefer to use detailed chemical kinetic models. There are two problems with this approach: 
1) Jet-A is a multi-component fuel and not a substance, and there are no universally accepted surrogate fuel models 
for Jet-A; 2) the computational costs associated with these models make them impractical when fine computational 
grids are used. Originally, a single-step, global chemistry model was used. This model was based on propane 
kinetics 21 , which are close to Jet-A’s reaction rates. The Jet-A fuel is modeled as single species (CnH 2 i) that 
represents a hydrocarbon mixture of decane, hexane and benzene in both gas-phase and liquid-phase solvers. The 
gas is treated as an ideal mixture with five-coefficient curve fits for C p (T) of each species and a CHEMKIN 
treatment of transport properties (species and ideal mixture rules). The single-step model allowed an easier start up 
in the solution process, by reducing the computational requirements during the ignition phase. Single-step models 
do not allow emissions calculations, only heat release. Because of this, a reduced twenty-step, fifteen species model 
was used based on the mechanism published in Ajmani-Kundu 22 . The reduced mechanism also describes the 
formation of Carbon Monoxide and Nitrogen Oxide. However, only one nitrogen-oxide species, NO, has been used 
in the reduced mechanism. NO in the reduced mechanism represents the whole family of nitrogen oxides including 
nitric oxide by Zeldovich 23 reactions, prompt NO reactions by Fenimore 24 , and nitrogen oxide formation through 
nitrous oxide. 

3. Liquid Phase Modeling 

The liquid Jet-A spray is evolved using a dilute spray Lagrangian solver (LSPRAY) 4 which neglects any 
drop-drop interactions. Turbulence effects follow the KIVA-II approach 25 of adding a turbulent fluctuation velocity 
to gas-phase velocity when calculating droplet drag and vaporization. The evaporation model includes solution of 
internal drop temperature distribution (thirteen point one dimensional mesh, finite -difference solution of a Hill 
vortex model) and a droplet regression rate employing three different correlations depending on droplet Reynolds 
number, Re drop (where Re drop is based on the relative speed V gas + V gas - V drop ). 

Jet-A fuel was represented as a single component, CnH 2 i. Fuel injection, atomization, and evaporation were 
handled in the following manner. Fuel injection in the model was handled by using Lefebvre’s 26 airblast correlations 
for a plain orifice. The correlation provided mean fuel droplets sizes of 16.3pm for the main fuel circuit and 13pm 
for the pilot. These particles are then injected, assuming a Rosin Rammler drop size distribution 27 , in the appropriate 
locations in the model, as jets in crossflow. These jets were presumed to have a solid cone angle of 7°, with each 
cone represented by eight circumferential droplets groups, with eight polar droplets groups at each circumferential 
location, stochasitally injected (with respect to polar and circumferential locations) every 200 spray solver time steps 
(40 microseconds). The main and pilot circuit drops were injected with velocities of 42 m/s and 13 m/s, respectively. 
Once in the flow, the particles were tracked 0.2 milliseconds before evaporation was switched on. 
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4. Computational Procedure 

The simulations are steady-state and solve the Favre-averaged transport equations for species, momentum and 
energy, as well as a two-equation k-s turbulence model. Staging was used in the solution process; cold-flow 
calculations and initial combustion calculations were performed using a single-step chemistry model with lagrangian 
spray until a steady state solution was obtained. The final stage of CFD calculations was performed by switching 
from the one-step chemistry model to the reduced chemistry model; this was done by changing the input chemistry - 
parameters of the code. It is important to note that no turbulence-chemistry interaction model was used for this case, 
so called “laminar chemistry”. So, averaged temperature and averaged mass fractions were used to compute the 
reaction rate. We believe this is appropriate as an engineering assumption, because for this particular case, turbulent 
kinetic energy (k) is below 20 m 2 /s 2 , with the peak on occurring at the tip of the fuel injector. 

Calculation of the source term due to chemical kinetics uses the explicit “reference species” approach described 
in the KIVA-II manual. In the current implementation, the explicit “reference species” integration of chemical 
kinetics source is performed using ten sub-iterations within each pseudo-time step of gas-phase Runge-Kutta 
integration in pseudo-time. 

The NCC computations for reacting and non-reacting flow were run in general until the flow residuals were 
reduced three orders of magnitude. The mass flow rates at the boundary conditions were also monitored as a 
convergence criterion. Dissipation (JST type) was set at 0.0 for second order dissipation (s 2 ) and 0.05 for fourth 
order dissipation (s 4 ). 28 The value of k 2 , the constant that scales the second order dissipation gradient switch, was set 
at 0.70. Setting the second order dissipation to zero is absolutely necessary to accurately resolve flow features, like 
jets. A CFL number of 1.0 was used. A cubic, non-linear k-epsilon model with a variable C mu coefficient was used. 
This model was selected because of the swirling flow. A dynamic wall function with pressure gradient effects was 
used to model near wall turbulent flow effects. 

Computations were performed on the SGI ICE computer “Pleiades” at NASA Ames. The “Pleiades” 
supercomputer with “Ivy Bridge” processors was preferred because its processor cores had considerably higher 
computational performance and because of its high speed, low-latency interconnect. This interconnect was important 
because the Lagrangian spray model created a load unbalance in compute nodes. The high speed interconnect 
seemed to mitigate the load imbalance. It takes approximately one week (wall clock) to complete a single PICS 
combustion case (starting from no initial solution) using 960 processors. 

IV. Results and Discussion 

Along with the emissions results obtained by gas sampling at the combustor exit, we present the CFD results for 
supersonic cruise and compare those with measured fuel patternation and Raman fuel and C0 2 concentrations 
relative to nitrogen. We also compare cases related to the ERA subsonic N+2 cycle for approach and elevated “idle” 
conditions, for fuel patternation via PLIF and PLS only. 


A. Supersonic Cruise 

Figures 6-9 show highlights of the CFD results at the N+3 supersonic cruise conditions run at NASA Glenn 
using fuel at ambient temperature. Later, we compare the simulated results with the laser-based measurements. The 
mean velocities and vorticity are shown in Fig. 6, while Fig. 7 shows the temperature distribution. Figs. 8 and 9 show 
the OH and fuel profiles, respectively. 

The PICS generates high radial velocity and vorticity within the swirler, which translates to a strong central 
recirculation zone immediately downstream from the swirler to help promote flame stability. The highest 
temperatures, from about 2200K - 2400K extend as far downstream as the recirculation zone and are bounded 
radially by the mid-level vorticity (with values on the order of 5500s" 1 ). 
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Figure 6. CFD results: normalized velocity and vorticity for N+ 3 supersonic cruise at center of combustor. Flow is 
from left to right, a) total velocity; b) axial velocity; c) radial velocity; d) vorticity. 


Figures 8 and 9 show the predictions for OH 
and fuel. We see high levels of OH — which show 
areas of reactivity and temperature — in the regions 
where fuel and air mix, so the regions of Jet- A (as 
represented by C11H21) coincide with high OH. 
Most of the OH occurs within the recirculation 
zone, but secondary and tertiary reactions persist 
downstream until reactions are complete. Fuel from 
the main stage projects downstream approximately 
to the axial location coinciding with the maximum 
reverse flow velocity and radially outside of the 
recirculation zone. This allows us to compare the 
fuel spray as measured using PLIF, PLS, and SRS 
to the CFD analysis. 

The optical measurements occurred beyond the 
PICS “dump plane” which is the injector’s 
downstream-most vertical end surface. Figures 10- 
15 show the results obtained at supersonic cruise 
conditions. The PLIF signal — representing fuel 



Figure 7. CFD results: Predicted contours of temperature 


within the Y=0 plane of the PICS combustor at the N+3 
supersonic cruise condition. Flow passes from left to right. 
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Figure 8. CFD results: Predicted contours of OH mole 
fraction within the Y=0 plane of the PICS combustor at 
the N+3 supersonic cruise condition. Flow passes from 
left to right 



U 26827E-06 
1.9307E-06 
1 3895E-06 
IE-06 


Figure 9. CFD results: Predicted contours of fuel mole 
fraction within the Y—0 plane of the PICS combustor at 
the N+3 supersonic cruise condition. Flow passes from 
left to right.. 


liquid and vapor — arises from the family of compounds having double ring aromatics (naphthalene — C i 0 H 8 ), as its 
base, such as methylnaphthalene. PLS signal arises from the Mie-type scatter presumedly of liquid phase fuel 
droplets on the order of 280 nm and larger. The SRS fuel signal is derived from the spectral region of the main C — 
H stretch, having a Raman shift near 3015 cm-1, near methane Vi and v 2 and ethylene Vi vibrational energies. The 
SRS signals for fuel and C0 2 are presented as mole fractions, normalized to molecular nitrogen. The PLIF and PLS 
signals are measured in counts of signal and presented on a linear scale. The signals are scaled separately to ease 
comparison. 

Figures 10, 11, and 14 show the fuel pattern for PLIF or SRS, each overlaid with CFD results at the 
corresponding spatial locations. The PLIF and PLS images were obtained on the same day, within minutes of each 
other; the SRS data were obtained on a different day. In each case, fuel signal emanates from a different “species”. 
Results from the optical techniques compare well with the numerical predictions, and show the fuel-air mixture is 
not uniform, which is to be expected because the fuel was injected at room temperature. There was not enough time 
for complete vaporization and mixing before exiting the main swirler. The amount of liquid compared to vapor can 
be seen qualitatively in Fig. 12, which shows the fuel PLIF signal (in red-yellow tones), comprised from both liquid 
and vapor signals; along with the PLS signal (blue tones), which is from the liquid fuel. Of the cases tested, this 
supersonic cruise point exhibited the highest signal from both PLIF and PLS. The two signals overlap almost totally 
to show that the fuel at this axial location is mostly in the liquid phase. Figure 13 shows the decay of fuel signal with 
axial distance from the dump plane for PLIF( solid red line) and PLS (dashed blue line). The mean signal was 
computed by calculating the average total signal at each axial location within the data blocks, then normalizing 
based on the respective maximum signal. The liquid signal decays at a faster rate than the total signal. We assume 
the fuel is effectively vaporized when the PLS signal reaches 5% of the maximum. Using this assumption, the fuel is 
evaporated an axial distance approximately 60% of the distance it takes the total remaining fuel to be consumed. 


Figures 11 and 14c show, respectively, the fuel PLIF and SRS signals, along with the CFD fuel in the y = 0 
plane, just downstream from the PICS dump plane. One may notice in Fig. 1 1 that the fuel PLIF signal contour is 
flat at the top; this marks the top of the optical field-of-view. All show similar results in that the fuel is projected 
outward and downstream from its origin. The slight displacement between measurements might be for a number of 
reasons. For example, the inlet conditions can vary slightly from one test day to the next. Setup of the two imaging 
techniques uses different equipment and necessarily different fiducials to designate a home reference position for the 
optics. SRS measurement results may also appear displaced because of the low spatial resolution of the 
measurement volume compared to the CFD or PLIF resolution; thus for example, the SRS fuel signal in Fig. 14c 
appears somewhat lower compared to the CFD. Also noted earlier, the test rig shifts depending on the inlet 
conditions, and some misalignment may occur based on rig growth. Although the CFD requires compromise in 
choices of models, mesh, and computational procedure, among others, it represents a physical ideal, while the 
research hardware, in particular, is likely to have some machining non -uniformities or imperfections and undergoes 
wear and tear for each test installation and as time progresses. 
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Figure 15 shows the results for carbon dioxide. The CFD mass fractions are shown as grayscale line contours, 
with the contours labeled. The SRS data are shown as a color scatter plot, with the symbols sized based on the 
concentrations relative to molecular nitrogen. From the CFD, we see that C0 2 concentration peaks in two regions: 
centrally, immediately downstream from the pilot; and in the annular region downstream from the main. The greater 
amount is near the pilot. For the Raman measurement, we found more C0 2 near the pilot, with decreasing amounts 
as one moves radially outward to the main region. The discrepancy of there being very little relative C0 2 
concentration near the main may be due to an increase in luminosity in the background. The increase in background 
luminosity for the top two rows of the detector is approximately thrice and twice that of the bottom three rows. The 
source of the increased luminosity might be flame chemiluminescence arising from the main flame that is just above 
the field of view, or from blackbody radiation emitted by the casting, which is also just above the field of view. This 



Distance from Center 

Figure 10. PICS Fuel pattern in an axial slice at the 
supersonic cruise condition: CFD prediction as line 
contour, with measured PLIF as color contour. Aft- 
looking-forward perspective. 


Axial Distance 

Figure 11. Comparison at supersonic cruise 
condition between predicted Jet- A (line contour) and 
measured fuel PLIF (color contour) distribution in 
the center vertical plane, at y = 0. Flow is from left to 
right. 



Figure 12. Supersonic cruise result: Comparison 
between total fuel signal represented by Fuel PLIF ( red- 
yellow contour ) and liquid fuel from PLS (blue contour ). 
Perspective view is aft-looking -forward. 



Figure 13. Supersonic crise result: Comparison 
showing the decay of liquid fuel (laser scatter, dashed, 
blue line ) to Fuel PLIF (liquid+vapor, solid red line) at 
supersonic cruise. Each plot is normalzed using its 
maximum signal. 
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increase in background luminosity decreased the signal to noise ratio to a level such that the signal was obscured. 
The plot in Fig. 15b is located 10mm downstream from Fig. 15a, and shows an increase in the amount of carbon 
dioxide including the area immediately downstream from the main annulus; the result is consistent with expectation 
and also with the thought that fuel interferences have some effect on the measurements in the region of the main 
injector. 
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Figure 14. Fuel fraction contours determined using CFD (grayscale, with lines) and SRS (color) at supersonic 
cruise condition, a) and b) are presented as axial slices, from the perspective of aft-looking-forward, with a) near 
the dump plane and b) 10 mm downstream. The 14c) image shows a center vertical slice with flow from left to right. 



Figure 15. Predicted and measured C0 2 signals within two axial slices, supersonic cruise condition, a) near the 
PICS dump plane, and b) 10 mm downstream. CFD results are presented as labeled, gray scale contours and the 
SRS results are scatter plots, wit size represented by color. The perspective is aft-looking-forward. 
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B. Subsonic N+2 cases 

NASA’s N+2 subsonic engine is expected to have a larger core than the supersonic N+3 engine, and higher 
operating pressures. Thus, it has higher air and fuel flows in its cycle. We provide a brief synopsis of results. 


1. Subsonic cruise 

Figure 16 shows an axial slice (aft-looking-forward 
perspective) near the dump plane, for the subsonic 
cruise condition (Table 2). We saw earlier that test 
conditions using fuel at ambient temperatures still 
exhibit a good deal of fuel in the liquid state after 
entering the combustor. However, for this condition, 
we determined using the method described above, as 
well as considering the total PLS and PLIF relative 
signals, that at most 10% of the total fuel signal was 
from the liquid. The higher inlet temperature and 
pressure, along with the commensurate higher flow 
rates and lower FAR contributed to higher 
vaporization and mixing, which is reflected in the 
more uniform fuel distribution compared to that shown 
in Fig. 12 for the supersonic cruise condition. This 
case also performed well for emissions, with 
EINOX < 5. 



Distance from Center 


Figure 16. Result for N+2 Subsonic Cruise. Fuel 
pattern near dump plane for subsonic cruise condtion as 
tested at NASA Glenn. The perspective is aft-looking- 
forward. 


2. Approach 

Figures 17 and 18 show results from two cases for approach conditions (Table 2), one with fuel staging between 
pilot and main (Approach- 1), the other using pilot only (Approach-2). The gas analysis results showed higher NOx 
levels for Approach-2 with EINOX 27% higher than for Approach- 1. Approach- 1 used approximately 8% more fuel 
than Approach-2, with all other inlet parameters equal. In this comparison, the local equivalence ratio was 
important. Figs. 17a and 18a show fuel patternation at the same axial position, using PLIF and PLS, and show that 
Approach-2 is less mixed than Approach- 1. We also noted for Approach- 1 that the PLS raw signal exceeded the 
PLIF signal; this is the only point tested for which this was the case, and indicates the majority of fuel was in the 
liquid phase. Figs. 18b and 19b show the corresponding C 2 * chemilumenescence, obtained by computing the mean 
from 200 instantaneous images, with exposure time for each image of 100 ps . Flow is left to right. These images 
give an indication of where much of the carbon-carbon breakup occurs. Reactions are occurring closer to the dome 
for Approach- 1 (implying reactions will be completed sooner), and without as much total intensity, than for 
Approach-2. These results indicate that if possible, fuel staging will be a better option for emissions during approach 
than running with pilot only. 
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a) 



axial position 


b) 

Figure 17. Results from Approach-1, a) Total fuel 
signal represented by Fuel PLIF ( red-yellow contour ) 
and liquid fuel from PLS (blue contour ), viewed aft- 
looking -forward. b) Mean C 2 * chemiluminescence 
signal within the vertical central plane (Y=0) obtained 
by computing the average of 200 individual images. 
Flow is left to right. 


Fuel PLIF 



axial position 


b) 

Figure 18. Results from Approach-2, a) Total fuel 
signal represented by Fuel PLIF ( red-yellow contour ) 
and liquid fuel from PLS (blue contour ), viewed aft- 
looking-forward. b)Mean Cf* chemiluminescence signal 
within the vertical central plane (Y=0) obtained by 
computing the average of 200 individual images. Flow 
is left to right. 


3. Idle 

The results at the elevated pressure idle condition (Table 2) are presented in a similar fashion as for the approach 
cases, in Figs. 19 (case Idle-1) and 20 (Idle-2). The idle cases use only the pilot and case Idle-1 has FAR 1.7 times 
case Idle-2, with both running fuel-rich. In this comparison, Idle-1 has less uniformity in fuel pattern and much more 
reactivity farther downstream, and resulted in higher EINOX levels. 
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Fuel PLIF 


Fuel PLIF 




a) 

Figure 19. Results from Idle-1. a) Total fuel signal 
represented by Fuel PLIF ( red-yellow contour ) and 
liquid fuel from PLS (blue contour ), viewed aft-looking- 
forward. b)Mean C 2 * chemiluminescence signal within 
the vertical central plane (Y=0) obtained by computing 
the average of 200 individual images. Flow is left to 
right. 



axial position 


b) 

Figure 20. Results from Idle-2. a) Total fuel signal 
represented by Fuel PLIF ( red-yellow contour ) and 
liquid fuel from PLS (blue contour ), viewed aft-looking- 
forward. b)Mean C 2 * chemiluminescence signal within 
the vertical central plane (Y=0) obtained by computing 
the average of 200 individual images. Flow is left to 
right. 


C. Emissions Measurements 

NOx emissions results from the NASA liquid-fuel tests are presented in Fig. 21, for PICS injector operation over 
a range of inlet conditions. For these emissions tests, combustor inlet pressure was varied from 150 psia to 255 psia 
(including the 175 psia supersonic cruise inlet pressure), and combustor inlet temperature was varied from 850 F to 
1004 F (about 80 F below the supersonic inlet temperature). The measured EINOx data for all of these conditions 
are plotted, and show only weak dependence on inlet pressure and temperature over this range. Combustion 
efficiency was greater than 99.95% for all data points shown. Based on the raw measured NOx data, the PICS 
injector appears to deliver about 5 EINOx at the supersonic cruise fuel -air-ratio, FAR/FAR S lto =1.1, when operated 
on liquid-fuel only (i.e. unheated fuel that is not vaporized) in the NASA rig, at the combustor inlet temperatures 
and pressures listed. 
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An improved estimate of supersonic cruise EINOx can be obtained from the correlation developed by NASA 
under the Advanced Subsonic Technology (AST) Program, for various fuel -lean aircraft combustor 
configurations 29,30 . In this correlation, EINOx scales with inlet pressure as P3 0 ' 595 (roughly square-root of P) and 
inlet temperature as exp(AT3/194 K), or exp(AT3/350 R) in English units. Using this correlation, EINOx data from 
175-180 psia are shown in Fig. 22a after scaling to the supersonic cruise combustor inlet temperature of 1087 F. The 
scaled data show an EINOx value of about 6 at FAR/FAR S lto = 1.1 for supersonic cruise. We note, however, that 
in the NASA rig the combustor residence time was 
significantly longer than expected in an actual 
engine application, as a result of the 12-inch 
distance between the PICS injector exit (the 
combustor dome) and the emissions -probe location. 

Thus, in an actual engine combustor we expect the 
liquid-fuel PICS NOx emissions to be well below 
the 6 EINOx “worst-case” scenario predicted here. 

Figure 22b plots EINOx values at 850 to 900° F 
inlet temperature, scaled to the 329 psia sea-level- 
takeoff pressure for this supersonic engine cycle. 

Thus, at 329 psia and 890°F inlet conditions, for 
sea-level takeoff, the liquid-fuel PICS configuration 
is predicted to deliver an EINOx value of about 3 (or 
less, given that the combustor residence time in the 
NASA rig is longer than expected in an aircraft 
combustor application). 

Finally, we note that the FAR values plotted in 
Figures 21 and 22 are based on a calculated airflow 
to the combustor. Two calculations are accounted 
for: First, because some air leaked around the 
combustor dome in the NASA rig (bypassing the 
primary combustion zone), and because this air 
leakage was not constant (it worsened during the 
course of testing), combustor airflow was calculated 
from the known effective area of the PICS injector 
and the measured pressure drop across the injector 
during testing. Thus, the leakage air (which did not 
participate in combustion) was ignored. Second, as 
in the UTRC testing 3 , the combustor airflow was set 
14% below the cycle specification (making FAR 
14% higher than the cycle specification for each 
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Figure 21. Measured NOx emissions from the single- 
injector PICS combustor testing at NASA. For these 
tests , unheated liquid Jet fuel was supplied to both the main 
and pilot circuits , and emissions samples were obtained at 
several different operating conditions, including variations 
in inlet temperature, pressure, and fuel/air ratio. Emissions 
are plotted against FAR in the combustor, as calculated 
from total fuel flow and combustor airflow, normalized by 
the FAR at sea-lev el- take off (SLTO) for this engine cycle. 
FAR/FAR slto = 1.10 at supersonic cruise, as indicated by 
the red lines. The blue lines indicate the FAR/FAR SLT0 = 
1.0 sea-level-takeoff condition. 

power condition) to account for an expected 14% 

liner cooling-air requirement for a full-scale annular PICS combustor in an engine. This essentially corrects for the 
fact that the NASA flametube rig was ceramic-lined and not air-cooled, and therefore has a lower total air 
requirement than an engine combustor. We note that reducing the liner cooling -air requirement below 14%, such as 
enabled by new materials or cooling technologies, will improve NOx performance (reduce emissions) by enabling 
increased process air for leaner combustion. The NOx data in all cases are consistent in being reported without 
correction to standard humidity, since high-altitude supersonic flight occurs in dry air, and the facilities at both 
NASA and UTRC supplied dried air to the combustion test rigs. 
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Figure 22. Predicted NOx emissions from the single-injector PICS combustor testing at NASA, based on 
scaling the measured results of Fig. 21 to cruise conditions (left-hand panel) and sea -level- takeoff conditions 
(right-hand-panel) for the supersonic cycle (Table 1) by applying the NASA correlation 29,30 for NOx emissions 
from lean-burn aircraft combustors. 


V. Conclusion 

We have presented computational and experimental analyses of the UTRC PICS injector, which has potential for 
application in an engine to meet the Pratt & Whitney N+3 supersonic cycle, or a subsonic N+2 engine cycle. NASA 
has defined technology targets for near term (called “N+l”, circa 2015), midterm (“N+2”, circa 2020) and far term 
(“N+3”, circa 2030) that specify realistic emissions and fuel efficiency goals for commercial aircraft. Experimental 
methods were also employed for select points of the subsonic cycle engine; cruise, approach, and idle with a slightly 
elevated inlet pressure. For the supersonic case, a PICS -based combustor was modeled using the National 
Combustion Code. 

Experiments at NASA employed gas analysis and a suite of laser -based measurement techniques to characterize 
the combustor flow downstream from the PICS dump plane. Gas analysis results showed that at both supersonic and 
subsonic cruise conditions, EINOx met the target goal of below 5. 

We saw that for the supersonic cruise point, use of ambient temperature fuel at the de-rated air inlet temperature, 
did not result in fully vaporized fuel before exiting the main fuel -air mixing passage. Results from the optical 
techniques compared well with the numerical predictions, and show the fuel-air mixture is not uniform, which is to 
be expected because the fuel was injected at room temperature. Analysis of fuel PLIF and PLS results indicate that 
the liquid fuel is evaporated an axial distance approximately 60% of the distance it takes the total remaining fuel to 
be consumed. Overall, the results show potential for meeting NASA’s supersonic cruise emissions goals with liquid 
fuel in all circuits. 

The subsonic cruise point performed well for emissions, with EINOX < 5. For this point we determined, that at 
most 10% of the total fuel signal was from the liquid. The higher inlet temperature and pressure, along with the 
commensurate higher flow rates and lower FAR contributed to a higher vaporization rate and better mixing, which is 
reflected in the more uniform fuel distribution compared to that for the supersonic cruise condition. 

For the approach points, PFIF and PFS results show that Approach-2 with fuel staging between pilot an main is 
better mixed than Approach- 1 which uses pilot only and that it produced lower NOx emissions. Analysis of the 
Approach- 1 case also indicated that the PFS raw signal exceeded the PFIF signal, indicating that the majority of fuel 
was in the liquid phase. 

The idle cases use only the pilot. Results indicated that idle case with the higher FAR exhibits less uniformity in 
fuel pattern and much more reactivity farther downstream. This also resulted in higher EINOX levels. 
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Using scaling methods previously developed at NASA, analysis of gas analysis results indicates that, at 329 psia 
and 890°F inlet conditions, for sea-level takeoff, the liquid-fuel PICS configuration is predicted to deliver an EINOx 
value of about 3. 
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